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Abstract 

We discuss an experiment in which two magnetic insulators that both show a field-induced 
magnetic-ordering transition are weakly coupled to one another and are placed into an ex- 
ternal magnetic field. If the respective magnetic states can be interpreted as phase coherent 
Bose-Einstein condensates of magnetic bosonic quasiparticles, one expects the occurrence 
of Josephson effects. For two identical systems, the resulting d.c. Josephson effect formally 
represents a constant quasiparticle Josephson current across the weak link, which turns out 
to be unobservable in an experiment. For magnetic insulators with different critical fields, a 
spontaneous alternating quasiparticle current develops with a leading oscillation frequency 
Wax. that is determined by the difference between the critical fields. As a result of the cou- 
pling, additional sidebands appear in the energy spectrum of the coupled device that would 
be absent without phase coherence. We discuss the primary conditions for such an effect to 
take place and conclude that its detection can be feasible for a proper choice of compounds 
with suitable and realistic material parameters. 

Keywords: Tunneling, Josephson effect. Quantized spin models. Macroscopic quantum 
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1. Introduction 

Quantum spin systems in solids have been a subject of intense research, both theoretically 
and experimentally. A number of such systems show magnetic-field induced phase transitions 
at zero temperature that have been interpreted as a Bose-Einstein condensation of magnetic 
bosonic quasiparticles [H-Isj. 

In insulating dimerized spin-1/2 systems, for example, one expects a condensation of 
triplet bosonic quasiparticles ("triplons") above a certain critical field -ffc 0, IH where the 
energy difference between the ground-state singlet and the lowest excited triplet states van- 
ishes due to the Zeeman splitting. In the language of magnetism, this condensation corre- 
sponds to a field-induced antiferromagnetic ordering that is associated with the appearance 
of a staggered magnetization [H, [s], 0|. However, the stability of such a condensate j§, [lo| 
and even the applicability of the BEC concept 111^ have been questioned. 
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Already known Bose-Einstein condensates (BEC) include superfluid helium 14, 
lute atomic gas clouds [iGj and pumped exciton-polariton and magnon condensates 
Several experimental proofs for the existence of a macroscopic quantum state have been 
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reported, e.g., from interference experiments [2lH23 



26| and Josephson effect s |27H32 



of these systems [ij, 115 
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from the observation of vortices [2^ 
, and even superfluid properties can be ascribed to most 
Such hallmarks of a true BEC state are all based on the 
existence of a macroscopic phase- coherent state, but none of these has ever been observed 
in any insulating quantum spin system. 

In this article we propose an experiment to probe the macroscopic phase coherence in 
insulating spin systems by observing the a.c. Josephson effect across a weak link. A success- 
ful detection of this effect would represent a direct proof for the existence of a macroscopic 
quantum state. We start with a discussion about the macroscopic wavefunction and the 
chemical potential of a triplon BEC in a magnetic insulator at zero temperature (sections |2] 
and |3]). We then consider the Josephson equations for a system of two weakly coupled triplon 
BEC's (section H]) and show how the a.c. Josephson effect can be identified in principle. In 
section \5\ we discuss several contraints for a system with realistic material parameters, and 
we finally conclude under realistic assumptions that such an experiment is indeed feasible. 



2. The macroscopic wavefunction of a triplon BEC 

The amplitude of the macroscopic wave function tp = ^/ne^'^ describing a Bose-Einstein 
condensate of triplons in a magnetic insulator at zero temperature is related to the density 
n, here defined as = N/N^ (with A"^ the total number of dimers) so that < n < 1. To first 
approximation and near the critical field HcAhis density is proportional to the longitudinal 
magnetic moment = g/jBnNd = g/J^sN jsl, l8|, where g is the Lande (^-factor and hb the 
Bohr magneton. As the equilibrium Mz is a well-defined quantity for a fixed value of the 
external magnetic field H in the thermodynamic limit, the triplon number A^ = M^/gfiB 
can be considered to be conserved. 

The phase (j) of the macroscopic wave function in magnetic insulators is closely related 
to (but not identical with) the angle ip of the transverse magnetic moments within the plane 
perpendicular to the main magnetic field H 0, E3, 38|. A plausibility argument for this fact 



is that the particle number A^ is canonically conjugate to the phase of the macroscopic 
wave function on the one hand, but also (via the proportionality of A^ to the longitudinal 
magnetic moment Mz and therefore to the associated component of the angular momentum 
Lz parallel to H) to the angle variable (p in the plane perpendicular to H on the other hand. 
In contrast to the time dependent global phase 0(t) of the macroscopic wavefunction, the 
angle (f between the transverse magnetic moments and the frame of reference given by the 
crystal lattice must be constant in time, because the Landau-Lifshitz-Gilbert equation does 
not allow for a perpetual spin precession of the staggered magnetization in the magnetic 
ground state. This apparent discrepancy can be easily resolved if one interprets the global 
phase of the macroscopic wave function as (f){t) = ip — ut (with u provided by the time- 
dependent Schrodinger equation), i.e., the in-plane angle ip of the staggered magnetization 
corresponds to an undetermined phase constant (pQ. This interpretation of </> and if ensures 
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that the expectation values of Lz and of the energy E per triplon are (L^)^ = h and 
{E)^ = hjj, respectively. 

In an axially symmetric spin system, (p must be able to take any arbitrary constant 



value between and 27r [37j, which is in certain contrast to experimental observations where 



if seems to be locked to a material-specific value [39|, |40|. We will address this issue of a 
violated axial symmetry later on in section 15.21 In the following we shall focus on the time 
evolution of the global phase (j){t) of the macroscopic wave function because its relation to 
the angle if does not play any further important role here. 

3. The chemical potential fi 

A decisive prerequisite for an a.c. Josephson effect to take place between two BEC is the 
presence of a steady non-zero difference A/x between the respective chemical potentials. In 
a condensate of magnetic quasiparticles in an insulating spin system, this chemical potential 
is usually taken as /i = gfxsfJ'oiH — He) [H, 0, ll]- As in all BEC of weakly interacting 
light bosons in the dilute limit n <^ 1, a non-zero n is related to n according to n = fi/vQ. 
The interaction constant vq describes the repulsive hard-core interaction between the bosonic 
quasiparticles, and is in magnetic insulators determined by the finite inter-dimer interactions 
which result in the formation of dispersive energy bands of the triplet excitations It is 
important to emphasize that the appearance of a non-zero value of fi = dE/dN (with E the 
total energy of the condensate) is common to all interacting BEC that are treated within 



the Gross-Pitaevskii formalism |41|, |42[ , and is not restricted to the magnetic systems under 
study here. 

In the following we aim to bring two of such condensates into contact in order to investi- 
gate possible analogues to the d.c. and a.c. Josephson effects. Therefore it is indispensable 
to examine whether or not the above definition of the chemical potential /i is actually ap- 
plicable to treat such an experiment correctly. If we consider, for example, two different 
magnetic insulators {a and fi) with different critical magnetic fields H^^a and i^c,/? that 
are placed into a common external magnetic field if, an apparently constant difference 
A/i = /io — = gfiofiB{Hc^j3 — Hc^a) is maintained by the external magnetic field (see Fig. 
[1]). This situation is formally equivalent to a device composed of two pieces of the same 
material with a single He-, but placed into two different magnetic fields differing by Aif . As 
we shall see later on, the necessary field gradient for the latter type of experiment to be 
successful is of the order of fio^H ^ 0.5 T along a length of ~ 0.5 nm. This is technically 
out of reach, and we therefore do not consider this scenario any further. 

This non-zero difference A/i is not simply a result of choosing different reference points 
of the energy scale for the two condensates. Firstly, the difference in the critical fields 
Hc,j {j = <y or P) originally stems from a difference in microscopic inter- and intra-dimer 
coupling constants that determine the individual energy gaps, i.e. the energy separation of 
the respective ground-state singlet and the triplet states in zero magnetic field. To close 
these energy gaps, different external magnetic fields Hcj must be applied beyond which the 
respective ground-state triplons condense. On an absolute energy scale (with the zero point 
chosen for a state with isolated spin-1/2 particles in vacuum and in zero magnetic field) the 
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Figure 1: Sketch of an experiment in which two magnetic insulators with different critical fields Hcj (j ~ 
a, P) are weakly coupled to one another through a magnetically inert spacer layer / and arc placed into an 
external magnetic field H > Hcj. The respective condensate densities tij and the chemical potentials jij 
vary approximately linearly with H. 



condensation occurs at the energy level of the respective singlet states, i.e., at Ej = —S/AJj 
per dimer, where Jj denote the respective intra-dimer coupling energies. Secondly and 
most importantly, the relevant quantity that enters the problem of treating the Josephson 
effect is the relative energy gain (loss) dE/dN upon creation (or annihilation) of one triplon 
quasiparticle to (or from) the condensate, from (or into) the singlet sea, which is actually 
UjVQj = jjij = gjjjBl^oiH — Hcj) per dimei0 jsl-lsl, 0|. Therefore, the jjj (and along with them 
the condensate densities rij = fij/voj) are indeed different in a common external magnetic 
field that exceeds both ifc,Q and Hc^p (see Fig. [1]), and A/i can then be considered as the 
analogue to an external voltage controlling the difference between the chemical potentials 
in a superconducting Josephson junction. 

4. Josephson effects 

4.I. The Josephson equations 

We now apply this concept to a system of two dimerized spin systems at zero temper- 
ature and with different critical fields i^c,/? > -f^c,Q beyond which the respective magnetic 
quasiparticles are supposed to condense. We initially assume a perfect axial symmetry of 
the two magnetic systems with respect to the direction of the external magnetic field, but 
we will discuss the case of violated axial symmetry later on in section 15.21 The boundary 



^This result for ^ can be illustrated by identifying E with the energy of the magnetized system in a field 
H,E^HoJI^^ M^{H') AH' with = giiBN and N = gi-iBl-i-a (H - He) Nd/vo, so that E{N) = N^va/2Nd 
and fi = dE/dN = uvq. 



4 



layers of the two systems are assumed to be weakly coupled to one another, and we place 
the device into a magnetic field H > Hcj (see Fig. [1]). For simplicity, we also assume 
■^^o.a ~ Vo^i3 = "^0 and Qa ^ Qp = 9, which is not essential for the main conclusion of our 
consideration, however. 

In an approach introduced by Feynman 43|] to explain the Josephson effects across 
a weak link between two superconductors one considers the macroscopic wave functions 

= y^/nje^'^^ on opposite sides of the junction, and treats the weak coupling between them 
according to ih-^ipa^fs = lia;p'^a;p + Kijjp-a, where K <^ is a phenomenological coupling 
constant. The resulting differential equations for the number of particles rij occupying the 
respective macroscopic quantum states and for the corresponding phase difference A0 = 
(pa — (pa become 



dria dni3 2K • /AA^ \ 

- - —y/n~npsm{A(p) , (la) 



dt dt h 

dA(p fia — fJ'13 K Ha — Up 
dt h h y/UaTlp 



COS (A0) . (lb) 



For Hc^j3 > Hc^ai the solution of Eq. [Tb] in the weak-coupling limit K <^ /i^ is, to first 
approximation, A0 ^ A^o + Wa.c.^ with Wa.c. ~ Afi/h and an undetermined constant phase 
difference A(pQ. The resulting variations in nj{t) are very small, i.e, by a factor K/Afi 
smaller than their time-averaged values. 

4-2. The d.c. and a.c. Josephson effects 

In the case of identical systems (A/x = and = np, see Fig. [2^) a constant A0(t) = 
A0O = (po,i3 — (po,a would represent a constant quasiparticle current dria/dt across the weak 
link (d.c. Josephson effect). As the constants 0oj correspond to the in-plane angles ipj of 
the respective transverse magnetic moments that will tend to align in the coupled device, 
we expect AcpQ = and therefore drij/dt = 0. 

In a device with two magnetic insulators with different critical fields, however, A/i 7^ 
and therefore cja.c. 7^ (Fig. [2t). As one would expect it from the analogy with super- 
conductors and superfiuids, the resulting variation dria/dt given by Eq. [Ta] represents an 
oscillation in time with the field-independent leading frequency Wa.c. (a.c. Josephson effect). 

The analogous situation in a superconducting Josephson device leads to the appearance 
of Josephson electrical currents. If connected to an external charge reservoir and fixing 
Afx to a constant value, the net charge-carrier density remains constant despite a nonzero 
dua/dt, but a measurable electrical current proportional to dn^/dt fiows from the charge 



reservoir and the attached leads through the whole device [43 1. 

With a junction composed of insulating spin systems it is the external magnetic field that 
maintains the difference between the chemical potentials to a certain value. This field fully 
penetrates the whole sample volume, and it therefore entirely "short-circuits" each of the 
two spin systems separately, thereby keeping the respective total numbers of quasiparticles 
Uj constant. As a consequence, there is no directional macroscopic quasiparticle current 
flowing within the two individual branches on both sides of the device. One can think 
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Figure 2: Sketch of the expected energy scheme (right panels) for a) two identical weakly coupled magnetic 
insulators, b) two uncoupled magnetic insulators with different critical fields, and c) two weakly coupled 
magnetic insulators with different critical fields (a.c. Josephson effect), all for H > Hcj and measured 
relative to the respective ground-state energies. In the case c), the ESR-activc mode Aq should split by a 
detectable amount hw^.c (sec text). Arrows in the left panels symbolize the quasiparticle currents related to 
dua/dt that can be ascribed to the spacer region (see text). 



of the magnetic field as a quasiparticle source (or sink) that replaces (or removes) those 
quasiparticles that are crossing the weak link. It is only by virtue of the phase difference 
between the two materials that a quasiparticle current can still be formally ascribed to the 
magnetically inert spacer region (see Fig. [2t), but it may be impossible to detect it directly 
in an experiment. 

4-3. Experimental manifestation of the a.c. Josephson effect 

To second-order approximation [K ^ with Hc,a 7^ -f^c,/? and therefore Ua 7^ nj^), Eq. 
[Tbl denotes a narrowband frequency modulation of the phases 

(j)a;l3it) ^ </)o,a;/3 - ^^t - ^ . sin(a;a.c.^ + A(/)o), (2) 

with a modulation frequency Wa.c. and a modulation index of the order of K/Afi. This 
modulation of (pjit) has profound consequences on the macroscopic wave functions tpj = 
^yfi-Q^'f'j describing the two condensates, because equally spaced sidebands should appear 
in the energy spectrum of the coupled device. This dynamic effect should manifest itself 
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in a corresponding splitting of all energies that are associated with transitions from the 
(condensed) ground state to (uncondensed) excited states (see Fig. [2t). Such a splitting is 
absent for two uncoupled magnetic insulators with different critical fields (Fig. |2]d) and, of 
course, also in a scenario where no macroscopic phase coherence is present at all. Therefore 
a successful experiment in a coupled system that can probe the occurrence of sidebands with 
separation /iWa.c. would represent a very strong experimental support for the existence of a 
state with macroscopic phase coherence. 

To be more specific, a high-resolution electron-spin resonance (ESR) measurement of 
the transitions between the ground state considered here, and the nearest excited triplet 
states (referred to as Aq in Ref. ji^] and Eq(Q) in Ref. jisf), or a one-magnon Raman 
experiment on the so-called -E'_(Q)-mode 46|, should reveal a characteristic splitting of the 
corresponding modes by fkjg,,c. which is related to the a.c. Josephson effect (see Fig. 
This is a central result of this work, and we shall discuss in the following whether or not the 
proposed a.c. Josephson effect can be observed in a device with realistic material parameters. 



5. Constraints in real systems 

5.1. The lifetime of phase coherence 

We may expect a successful experiment only if the separation /jo^a.c. between associated 
sidebands in the energy spectrum is comparable to or larger than their linewidth which can 



be related to the inverse lifetime r^^^ of phase coherence and to the inverse lifetime r^^^ of the 

to the best 



> 



_ max[Tp^ 



T, 



qp 



The lifetime r^c has. 



magnetic quasiparticles involved, 
of our knowledge, not yet been estimated for triplon EEC's, but we may initially compare 
it to Tqp, i.e., that of the S = 1 quasiparticles. Scattering processes on magnons, phonons 
or impurities can be important limiting factors, but there is no fundamental principle that 
would restrict it to inaccessibly short time scales at low enough temperatures. Corresponding 
values for the quasiparticle lifetime Tqp ranging from > 10^^^ s up to ~ 5 x 10"^^ s can be 
inferred from precise inelastic neutron-scattering measurements of the magnon linewidth in 
the spin-dimerized compounds TlCuCla 



471 and (C4Hi2N2)Cu2Cl6 [48 



respectively, taken 

at low temperatures for quasiparticles that are not par t of the condensate, and from the 
ESR linewidth of the Aq mode observed in TlCuCla |44l |. 

It is conceivable that the actual lifetime Tpc of phase coherence in the BEC is even 
considerably longer than that of its constituting particles r^p, e.g., thanks to the formation 
of energetic barriers of topologic origin that may lead to a strong suppression of dissipation. 



i.e., to the occurrence of true spin superfluidity [37|. It is worth mentioning here that the 



a.c. Josephson oscillations that have been seen in an experiment on an exciton-polariton 
Josephson junction were clearly observable although the oscillation frequency (wa.c. ~ 4 x 



10 



was very near the inverse lifetime of the respective bosonic quasiparticles 



r. 



3 X 10" s 



32 



qp 



but still much larger than the inverse of the surprisingly long measured 
-1 



phase-coherence time (r Ri7xl0s j49 



I.e., r, 



>> Tqp). To be on the safe side we 



conservatively choose Wa.c.^ 10^^ s ^ > t^^ for our proposed experiment, which corresponds 

to flo{Hc,f3 - Hc,a) ~ 0.5 T. 
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5.2. Violated axial symmetry 

So far we have considered only perfectly axially symmetric materials. Any violation of 
axial symmetry of the magnetic exchange interaction leads to the occurrence of an anisotropy 



thereby lifting the Goldstone linearity for /c — )■ 0, 
H, [ij and leading to well-defined magnetic structures 



gap A in the energy spectrum E{k) [54 [, 
fixing the phases (poj to constant values 

with in-plane angles ipj that are locked to the crystal lattice 
estimated to A{H) = 8'yjj,{H), where 7 > is a measure for the exchange- interaction 
anisotropy in the plane perpendicular to H 5J|. For experimental energies larger than 



39l . l40| . This gap has been 



A, however, a qua^i-linear doldstone mode is^^overed 0, the effects of axial anisotropy 
are expected to be smeared out on short enough time scales, and we can ascribe an upper 
limit T/\ = h/ A to the lifetime of phase coherence so that we must additionally require 



5.3. The weak link 

In a real experiment, A/i = g^Q^B^Hc^p — i^c,a) and the coupling constant K are deter- 
mined by the junction properties. The dimensions of the weak link should not significantly 
exceed the healing length of the condensate (the analogue to the coherence length in su- 
perconductors). For weakly interacting Bose gases, this healing length is ^ = h/ ^/Inv^m* 
(where m* is the effective mass ges as the temperature 

is increased towards the Bose-Einstein condensation temperature where n vanishes. 

The simplest version of such a weak link may consist of a slit that is narrow enough to 
allow for a small but finite magnetic coupling and therefore for a tunnelling of magnetic 
quasiparticles. Such a junction might also be manufactured by the layered growth of a 
magnetic insulator (a, with lic,a) containing a certain species of magnetic ions forming spin- 
1/2 dimers, followed by one or a few layers of an isostructural compound with nonmagnetic 
ions, and by further adjacent layers of another isostructural magnetic insulator (/3, with 

5.^. A realistic toy model: (Ba,Sr)-^Cr20^ 

To obtain numeric values, we consider a model device composed of Ba3Cr208 (a) with 
l-i'oHca ~ 12T [i^], separated by nonmagnetic isostructural Ba3V208 51| from a compound 



Ba3_a;Sra.Cr208 (/9) [52p with x chosen in such a way to achieve a larger p,nHr r, but still 
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much smaller than /io-f^c2 = 30T of fully stoichiometric Sr3Cr208 with x = 1 [52] • If we 
achieve a fioHc,i3 ~ 12.5 T, the expected characteristic frequency becomes Wa.c. 
With an external magnetic field fioH ^ 13 T, and taking Vo/ks ~ 8.7K with g 



2 53 



we 



have ria ~ 0.15 and n/3 ~ 0.08. The choice of a small H — H^j (i.e., fij <^ vq) is primarily 
necessary to fulfill the dilute-limit condition rij <C 1, but it may not be crucial for observing 
the essence of the predicted effect. With m* ^ 1.5 x 10"^"^ kg and a mean distance d ^ 0.6 nm 



between the dimers [40|, |53|, the resulting values for the healing lengths at T = are small 



^We have already verified the existence of the solid-solution scries Ba3_:cSr:cCr208 for polycrystalline 
samples, but we have no further information about their magnetic properties. 
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i^a ~ 0.4 nm and ~ 0.6 nm ^ d), but a very narrow slit or a weak link composed of one or 
only a few unit cells of the nonmagnetic compound could be adequate. If still larger values 
of C,j are required, rij can be diminished by further reducing A/i and fip. Alternatively, 
the experiment can be performed near the condensation temperature as it has been done 
in superfluid ^He to make the a.c. Josephson effect observable j23], possibly at the cost of 
reducing the lifetime Tgp of the triplon quasiparticles. 

With the material parameters used for the experiment described above and cJa.c.~ 
10^^ s~^, 7 should be below ^ 4/ieV. The exchange-interaction anisotropy has been esti- 
mated to 7 f« 16 - 30 /ieV in TlCuCk 0, Issj, but to only 7 ^ 1 /leV < 4/ieV in BaCuSiaOg 



56| and particularly in Ba3Cr208 j40|, so that the proposed device would fulfil all of the 



requirements listed above. 



6. Concluding remarks 



To conclude these numeric estimates we want to state that other material properties, such 
as the presence of a strong Dzyaloshinski-Moriya interaction that we have not considered 



here, may still inhibit the formation of a phase coherent condensate [5J| and therefore make 
the observation of the a.c. Josephson effect impossible. However, we classify such factors 
as not intrinsic to the problem, and a proper choice of compounds with suitable material 
parameters should make the generation and the observation of the a.c. Josephson effect in 
magnetic insulators feasible in principle. Although we have chosen dimerized spin systems 
operating near their lower critical fields as a model for the present consideration, all the 
above arguments should also hold near the respective saturation fields as a consequence of the 
particle-hole symmetry of the problem 0] , and for all other types of insulating spin systems 
that are supposed to show a field-induced BEC of magnetic quasiparticles. It is conceivable 
that Josephson- like phenomena can even naturally occur in certain quantum magnets, e.g., in 
a system showing an intrinsic modulation of the boson density (and therefore of the chemical 
potential) due to inequivalent planes hosting the triplon condensate. Such a situation has 
been reported for BaCuSi208, in which the average boson density strongly varies along the 
c-axis js^, Esj . 

We finally mention here an interesting analogy between the case of an ESR measurement 
on weakly coupled magnetic insulators as discussed above, and corresponding microwave 
stimulated experiments on ferromagnetic films which are separated by normal-metal spacers 
(59I 1 . In both cases, there exists a dynamic coupling between two different materials through 
a magnetically inert spacer layer, thereby altering the dynamics in both systems as compared 
to the uncoupled situation. 

In summary, we suggest that two weakly linked magnetic insulators with different critical 
fields Hcj that are placed into a suitably chosen external magnetic field should show an 
altered energy spectrum as compared to the uncoupled limit, with additional sidebands 
that are separated by /kJa.c. = fi'/^_ByLto(i^c,/3 — -f^c.a)- These sidebands are a manifestation 
of the a.c. Josephson effect and can be tested, for example, in a high-resolution ESR 
experiment. Assuming realistic material parameters and considering several constraints, we 
conclude that this effect can indeed take place in a real system. A corresponding successful 
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experiment would represent a very strong experimental support for the existence of a state 
with macroscopic phase coherence. 
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